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Abstract 

In the present paper, we start from the canonical theory of loop quantum gravity and the master constraint programme. The 
physical inner product is expressed by using the group averaging technique for a single self-adjoint master constraint operator. 
By the standard technique of skeletonization and the coherent state path-integral, we derive a path-integral formula from the 
group averaging for the master constraint operator. Our derivation in the present paper suggests there exists a direct link 
connecting the canonical Loop quantum gravity with a path-integral quantization or a spin-foam model of General Relativity 
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1 Introduction 



Loop quantum gravity (LQG) is a mathematically rigorous quantization of general relativity (GR) that preserves background inde- 
pendence — for reviews, see |[Tl|2l[3). It is inspired by the formulation of GR as a canonical dynamical theory of connection pi]. In 
this formulation, the canonical coordinate on the phase space of GR is the su(2)-connections A'^ and the densitized triads E". The 
total Hamiltonian of GR is a linear combination of the Gauss constraint, the spatial diffeomorphism constraint and the Hamiltonian 
constraint. Thus the dynamics of GR are essentially the gauge transformations generated by the constraints. 

At the kinematical level of LQG, we smear the su(2)-connections A'^ along paths e to define the holonomies A(e), and smear the 
densitized triads on 2-surfaces to define the fluxes £'(5'). We consider an arbitrary finite piecewise analytic graph y embedded in 
the spatial manifold Z, which is a collection of edges (piecewise analytic paths) and their endpoints, called vertices. We denote £(7) 
and y(y) the collections of the edges and vertices of y respectively. A cylindrical function fy is a continuous function depending 
on the holonomies along the edges in the graph y. The collection of all the cylindrical functions Cyl with a certain sup-norm 
turns out to be an Abelian C*-algebra. The spectrum of this C*-algebra is a compact Hausdorff space f5l, which is denoted by 
^ and called the quantum configuration space. The C*-algebra of cylindrical functions is the algebra of continuous functions on 
the quantum configuration space Moreover it turns out that the fluxes E(S ) can be considered as the algebraic vector fields 
on y[. The collections of the cylindrical functions and the flux vector fields generates the so called, holonomy-flux *-algebra 31. 
The uniqueness theorem guarantees that there exists a unique state (positive linear functional) on this holonomy-flux *-algebra 91 
such that it is invariant under both SU{2) gauge transformations and spatial diffeomorphism transformations |6|. This invariant 
state gives a unique representation Hilbert space "Ha:,;]!] which is called the kinematical Hilbert space of LQG. This kinematical 
Hilbert space is isomorphic to the space of square-integrable function on yi, L^{J[, djiAi), where the measure hal on J[. is called 
Ashtekar-Lewandowski measure. This is the reason why we call J[ the quantum configuration space. Moreover, the space of 
cylindrical functions Cyl is dense in the kinematical Hilbert space 'Ka:,„. A special class of cylindrical functions, which is called 
the spin-network function T s=(y,j,m,n), forms an orthonormal basis in ^Kin- Since the spin-network function Ts=(yj.m,n) is labeled by 
a continuous parameter y, the kinematical Hilbert space 'HKin is non-separable. Furthermore, the Gauss constraint can be quantized 
and solved, which results in the Hilbert space of S U(2) gauge invariant states. The finite spatial diflFeomorphisms can be represented 
as unitary operators on 'HKin- Thus by using the group averaging technique, we obtain the Hilbert space of diffeomorphism states 
'Hoiff- All above is the beautiful kinematical framework of LQG, see [1 1 for all the detailed constructions. 

The dynamics of LQG is determined by the quantization of Hamiltonian constraint of GR. There are several different approaches 
to implement the Hamiltonian constraint at the quantum level. First of all, a quantum Hamiltonian constraint operator H{N) was 
first defined in Q. This operator is a densely defined closed operator on the kinematical Hilbert space 'HKin and featured by its 
graph-changing operation on any cylindrical function fy. On the other hand, given a number of first-class constraints Cj, we can 
define a, so called, master constraint by 



where K is a positive definite matrix. Correspondingly, an operator M is defined and usually, it is a positive self-adjoint operator. 
There are two versions of the master constraint operator in LQG. One is a positive self-adjoint operator defined on the Hilbert 
space 'Hoiff of diffeomorphism invariant states, with the graph-changing operations |8 |. The other version is non-graph-changing, 
positive, self-adjoint operator, defined on the kinematical Hilbert space l-lKm- This non-graph changing master constraint operator 
is also adapted in the framework of algebraic quantum gravity (AQG) 111, which define the quantum theory of GR on a single cubic 
algebraic graph y with an infinite number of edges and vertices. Its kinematical Hilbert is a infinite tensor product extension ifTOl of 
the LQG kinematical Hilbert space 'Hxin- It is remarkable that the semiclassical limit of the non-graph changing master constraint 
is correctly reproduced by the complexifier coherent states ifTTl [121 [T3J . With the Hamiltonian constraint operator or the master 
constraint operator, a Quantum Einstein Equation is well-defined as the quantum constraint equation H{NyV = or M*!' = cH in 
LQG, whose space of solutions with a chosen physical inner product is called physical Hilbert space ^Hphys- 

'in a way similar to that tlie states in Fock space are excitations from tlie unique (Poincara invariant) vacuum state. 
^See section[2]for the definition of M. 
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In order to rigorously define the physical Hilbert space from a self-adjoint master constraint operator, we follow the general 
procedure of direct integral decomposition (DID) lfT4l , which rigorously specify the physical inner product in general. This procedure 
is also tested in various physical models |fT5l and gives correct results. Another approach to obtain the physical inner product is the 
group averaging technique for rigging inner product (T6\ [T7] [TS] . In order to employ the group averaging technique for a given 
first-class constrained system, the constraint algebra formed by the first-class constraint should have a Lie algebraic structure so 
that the gauge transformations generated by the constraints form a group. In the case of GR, it is well-known that the constraint 
algebra formed by the diffeomorphism constraint and Hamiltonian constraint is not a Lie algebra, because of the presence of structure 
function. Thus the group averaging doesn't work directly for the diffeomorphism-Hamiltonian constraint algebra. However if we 
quantize the constraint system using a single master constraint, the constraint algebra is trivial thus is a trivial Lie algebra (since 
there is only one constraint). Therefore the group averaging can also be carried out for the master constraint, which has been akeady 
pointed out in fS) long time ago. In addition, we have shown that the group averaging for the master constraint gives the consisent 
result as it is given by DID procedure under some technical assumptions ifTSl . which are fulfilled by all the physical models tested 
in lfT4l . So if the master constraint operator in LQG satisfies those assumptions, the group averaging technique is the correct way to 
compute the physical inner product in LQG. 

The general consideration in the last paragraph is a preparation for linking canonical LQG with a path-integral formula. The path- 
integral formulation is of interest in the quantization of general relativity (GR), a theory where space-time covariance plays a key 
role. Currently the spin-foam model IfSOl can be thought of as a path-integral framework for loop quantum gravity, directly motivated 
from the ideas of path-integral adapted to reparametrization-invariant theories. The current spin foam approach is independent from 
the dynamical theory of canonical LQG [7 1 because the dynamics of canonical LQG is rather complicated, it uses an apparently much 
simpler starting point: Namely, in the Plebanski formulation 1.22 J , GR can be considered as a constrained BF theory and treating the 
so called simplicity constraints as a perturbation of BF theory, one can make use of the powerful toolbox that come with topological 
QFT's on . It is an unanswered question, however, and one of the most active research topics momentariljH, how canonical LQG 
and spin foams fit together It is one the aims of this paper to make a contribution towards answering this question. 

There are some early attempts in looking for a link between canonical LQG and a path-integral formulation of LQG. In the 
seminal paper [23], the authors presented a heuristic method in solving the (Euchdean) Hamiltonian constraint operator in LQG, i.e. 
taking any diffeomorphism invariant state if/ e "J^oiff and map to 



where H' is the dual operator on Cyl*. Eq.( |L2t have a shape similar as a group averaging. However such a method is formal 
because (1.) The exponential is not well-defined because H'(x) is not a self-adjoint operator; (2.) the product the operator H'{x) 
doesn't preserve the 'Hoiff- If we use W{x) and operating on kinematical state if/ e 'HKin, the product flxei, is not well-defined since 
for different jc € S, W(x) doesn't commute. However if we proceed formally, we obtain the physical inner product for given two 
diffeomorphism invariant spin-network states 



The physical inner product so defined has a spin-foam-like structure since the Hamiltonian constraint operator W(N) is graph- 
changing, i.e. add several arcs in the neighborhood of each vertex whose valence more than 3. But as we have seen, although this 
physical inner product suggests the spin-foam as the path-integral quantization of GR, it is unfortunately ill-defined, essentially by 
the fact that the constraint algebra formed by the Hamiltonian constraint is not a Lie-algebra. 

The present paper applies the general considerations in IfTTlfTSlI to the case of GR. Here we improve the situation in Eq.( ll.2l ) by 
using a single self-adjoint master constraint M. As it is shown in flSl, the group averaging using the master constraint operator is 
well-defined and correctly specifies the physical inner product and physical Hilbert space under some technical assumption^ In this 
paper, we will show that this group averaging physical inner product can be computed practically (under some technical assumptions) 

^Here we are refeiTing the spin-foam model for 4-dimensional gravity. 

'^More precisely, the group averaging using the master constraint operator correctly specifies the absolutely continuous sector of the physical Hilbert space under 
some technical assumptions. 
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by using the available semiclassical technique lfT2l . It turns out that the result of the computation lead to a path-integral formula of 
LQG, which is a discrete analog of the naive Ansatz: 



(1.4) 



up to a local measure factor when the matrix K in Eg. (II. lb is phase space dependent. 
The present paper is organized as the follows: 

In section |2] we briefly review the framework of algebraic quantum gravity and the master constraint programme, the definition 
of group averaging technique and its consistency with direct integral decomposition. We also briefly review the semiclassical tools 
in AQG, which will be employed in derive the path-integral. 

In section |3] we start from the definition of the group averaging physical inner product using the master constraint. By using the 
standard technique of skeletonization and coherent state path-integral, we derive a path-integral formula from the discrete setting. 
With some certain assumptions and approximations, we obtain a path-integral formula in analogy with Eq. (ll.4b up to a local measure 
factor. 

In section m we summarize and conclude. 



2 Algebraic quantum gravity and master constraint programme 
2.1 The framework of algebraic quantum gravity 

Algebraic quantum gravity (AQG) lID is a modified approach of loop quantum gravity. But in contrast to loop quantum gravity, the 
quantum kinematics of AQG is determined by an abstract *-algebra generated by a countable set of elementary operators labeled 
by a single algebraic graph y with countably infinite number of edge^ while in loop quantum gravity the elementary operators are 
labeled by a collection of embedded graphs with a finite number of edges. Thus one can expect that in AQG, we don't consider the 
information of the topological and differential structure of the manifold in all the quantization procedure except the semi-classical 
analysis. When we consider the semiclassical limit of AQG, we should specify an embedding X : y — > Z, which makes the contact 
of the abstract operators with the physical fields on the spatial manifold E. 

The quantization in AQG bases on a single algebraic graph, which only contains the information of the number of vertices and 
their oriented valence. 

Definition 2.1. An oriented algebraic graph y is an abstract graph specified by its adjacency matrix y, which is an N X N matrix. 
One of its entries yu stand for the number of edges that start at vertex I and end at vertex J. The valence of the vertex I is given by 
V/ = Yijiyij + yji)- We also use V{y) and E{y) to denote the sets of vertices and edges respectively. 

In our quantization procedure, we fix a specific cubic algebraic graph with a countably infinite number of edges = N and the 
valence of each vertex vj - 2x dim(2;). 

Given the algebraic graph y, we define a quantum *-algebra by associating with each edge e an element A(e) of a compact, 
connected, semisimple Lie group G and an element Ej(e) take value in its Lie algebra g. These elements for all e e E(y) are subject 
to the canonical commutation relations 

[A(e),A(e')]^0, 

[Ej(e),A{e')] = ihQ\yTj/2A(e), 

[Ej(e),Eie')] = -ihQ\,,,,fjk,E,(e'), (2.1) 

and *-relations 

A(er = [Mey'V, E/eT = Ej(e), (2.2) 

'As it will explained later in this subsection, the Hilbert space on the single graph y contain rich enough information, more precisely the Hilbert space is 
non-separable, since the graph y is infinite. 
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where Q stands for the couphng constant (Q^ - Kin the case of GR), tj is the generators in the Lie algebra g and fjki is the structure 
constant of g. We denote the abstract quantum *-algebra generated by above elements and relations by 91. 
A natural representation of 91 is the infinite tensor product Hilbert space IfTOl 



(2.3) 



where 'T/f = L^(G,t/yUH), whose element is denoted by®/ = ®efe- Two elements ®/ and®/- in '7/*' are said to be strongly equivalent if 
Tie I < fe^fe >H. ~ ^ I couvcrges. We denote by [/] the strongly equivalence class containing ® y . It turns out that two elements in "K® 
are orthogonal if they lie in different strongly equivalence classes. Hence the infinite tensor Hilbert space TY® can be decomposed as 
a direct sum of the Hilbert subspaces (sectors) "TY^j which are the closure of strongly equivalence classes [/]. Furthermore, although 
each sector 'K^j is separable and has a natural Fock space structure, the whole Hilbert space 'K® is non-separable since there are 
uncountably infinite number of strongly equivalence classes in it. Our basic elements in the quantum algebra are represented on 
in an obvious way 



A(e) 
Ejie) 



[A{e)fe] ® [^e'^efe'l 



(2.4) 



where is the right-invariant vector field on the Lie group G. As one might have expected, all these operators are densely defined 
and E j(e) is essentially self-adjoint. Given a vertex v e V(y), the volume operator can be constructed by using the operators we just 
defined 



y fine2nf3=i' 



(2.5) 



where the values of e,,(ei ,62,^3) should be assigned once for all for each vertex v e Vij). When we embed the algebraic graph y into 
some manifold, the embedding should be consistent with the assigned values of e, (ei , 62, e^). 

Then we discuss the quantum dynamics. The half densitized constraints can be quantized to be composite operators as we list 
below m. 



Gauss constraint 



^ Ejie)- 2 Oj,[Aie)]Etie) 

b(e)=v f(e)=v 



(2.6) 



• Spatial diffeomorphism constraint 

- y 



e,.(ei,e2,e3) 



EW) ^ |L(v, 61,62)1 



A(J3)-A(J3r 



A(e3) 



A(e3)- 



(2.7) 



• Euclidean Hamiltonian constraint (up to an overall factor) 

e,,(ei,e2,e3) 



A(I3)-A(I3)- 



A(e3) 



(2.8) 



• Lorentzian Hamiltonian constraint (up to an overall factor) 



f{v) Yj ^v{ex,e2,ei)tx(A{e,) A(e,r\[H^^\v\ A{e2) A{e2r\[H^^\v\ 



A{e,) 



Aie^r 



H 



(1/2) 



(v) = ff'l^\v) + t{v)- 



(2.9) 



where the matrix Ojk[g\ is the adjoint representation of § e G on the Lie algebra g, V := 2,, V,,, i/^'' '■- Tjv ^^V^) and 

Qt^ - Yj e,,(ei,e2,e3)tr(A(ei)[A(e,)-\yW] A(e2)[A(e2)-\y<:'] A^^^^ (2.10) 
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^(v') denotes the binomial coefficient which comes from the averaging with respect to the triples of edges meeting at given vertex 
V. L(v, 61,62) denotes the set of minimal loops starting at v along ei and ending at v along gj'. And a loop /3 e L(v,6i,62) is 
said to be minimal provided that there is no other loop within y satisfying the same restrictions with fewer edges traversed. Note 
that since we only have a single cubic algebraic graph, the diffeomorphism constraint can only be implemented by defining the 
operators corresponding to diffeomorphism generators because a finite diffeomorphism transformation is not meaningful in our 
algebraic treatment unless the algebraic graph is embedded in a manifold. As a result, the (extended) master constraint operator can 
be expressed as a quadratic combination: 



It is trivial to see that all the above operators are non-graph-changingj and embedding independent because we have only worked on 
a single algebraic graph so far. It is obvious that the master constraint operator M is a positive and symmetric operator, thus has the 
Friedrichs self -joint extension ||26ll . We take this self-adjoint extension of M and consider M as a self-adjoint operator 

Note that the master constraint operator M is not densely defined on the whole infinite tensor product Hilbert space 'Hku, ■ In order 
to make M densely defined we have to restrict ourselves in a subspace such that M is densely defined on the subspace. On the other 
hand, it is remarkable that the master constraint operator M preserves all the strong equivalence class sectors, which are separable 
Hilbert spaces although the whole 'Hkui is non-separable. Because of this property, we can use the direct integral decomposition 
(DID) separately in each strong equivalence class sector to define the physical Hilbert space for each sector (see the next section). 
Then the total physical Hilbert space is a direct sum of the sectorial physical Hilbert spaces. Moreover, the strong equivalence class 
sector "Hal where the vacuume state at - I lives is especially interesting. A generic states in the sector "Hal only has excitations on 
a finite number of edges, which shows the similarity with the kinematical Hilbert space in LQG. 

2.2 Physical inner product and group averaging technique 

Given the self-adjoint master constraint operator M in Eq.( l2.1 11 1, we can formally define the quantum master constraint equation by 



The space of solutions for this equation combined with a certain physical inner product is called the physical Hilbert space 'Hphys- 
However, firstly the equation Eq.( l2.12b is only formally written because zero is probably contained in the continuous spectrum of 
the master constraint operator, so that the solution state *P may be not living in the kinematical Hilbert space anymore. Secondly it is 
not clear from Eq.( |2.12| i about which physical inner product should be chosen on the space of solutions, unless zero is a pure point 
spectrum and the space of solutions is a subspace of 'HKin- In order to rigorously define the physical Hilbert space "Hphys, we should 
in principle employ the direct integral decomposition (DID) [ l4| for the master constraint operator M. Note that DID requires that 
the Hilbert space should be separable. Fortunately, since the master constraint operator M is self-adjoint and preserves all the strong 
equivalence class sector in the infinite tensor product Hilbert space, the physical Hilbert space "Hpiivs is well-defined in principle for 
each strong equivalence class sector Finally we should make a direct sum for all the resulting sectorial physics Hilbert spaces. 

In lfT4l . the programme of direct integral decomposition (DID) is introduced in order to rigorously define the physical Hilbert 
space for a general constraint system. Such a programme proceeds as the follows: 

1. Given a separable kinematical Hilbert space 'Hku, and a self-adjoint master constraint operator M = K'^c\Cj, first of all we 
have to split the kinematical Hilbert space into three mutually orthogonal sectors 'HKin - 'H'''' ®'H"'^ Ql-C^^ with respect to the 
three different possible spectral types of the master constraint operator M. 

2. We make the direct integral decomposition of each li*, * — pp, ac, cj^l with respect to the spectrum of the master constraint 
operator M restricted in each sector, i.e. 



The non-graph-changing operator means that the operator doesn't change the underHne infinite cubic graph y, i.e. it doesn't add any edge in addition to the edges 
in E(y). 

^It labels respectively pure point, absolutely continuous, or continuous singular spectrum of the self-adjoint operator. These names come from the properties of 
the spectral measures on the coiTcsponding sectors. See 1141 for detailed definition. 



M := J] [G(.'/2>(v)^Gy^'V) + d\"^\v)^D^I'^\v) + H^'l^\v)^ H'^'I^\v)]. 



(2.11) 



M*!* = 



(2.12) 
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3. Finally we define the physical Hilbert space to be a direct sum of three fiber Hilbert spaces atA-O with respect to three kinds 

^,)=0 ^ '".1=0 ^ '".1=0 



of spectral types, i.e. Uphys = "HZ,-, © 'H"' ® 'HI 



Note that in step 2, we have assumed that all the ambiguities outlined in lfT4ll have been solved by considering some physical 
criterions e.g. the physical Hilbert space should admit sufficiently many semiclassical states, and it should irreducibly represent the 
algebra of Dirac observables as an algebra of self-adjoint operators. With this assumption, the procedure of DID programme gives 
a precise definition of the physical Hilbert space for a general constraint system. And in many models simpler than GR, such a 
programme gives satisfactory physical Hilbert space 1 15 1. 

However, if we want to practically obtain the physical Hilbert space of AQG and the detailed knowledges about the structure 
of this physical Hilbert space, then DID is not a suitable procedure. The reason is the following: The whole procedure of DID 
depends on a good knowledge of the spectral structure for the master constraint operator While for the case of LQG or AQG with a 
complicated master constraint operator M, the spectrum of M is largely obscure so that the DID programme is too hard to proceed 
practically even the first step. Therefore for our practical purpose in exploring the AQG physical Hilbert space, we have to employ a 
technique such that the final structure of physical Hilbert space l-lphys - '^jfo ® '^"=0 ® '^1=0 ^'^^ some sectors in it) is written down 
at least formally without much of the knowledges for the spectrum of the master constraint operator. Fortunately we have a single 
constraint in the quantum theory, whose generated gauge transformations for a one-parameter group. Therefore we can employ an 
alternative, (modified) group averaging technique to obtain the physical inner product IfTSl : 

Definition 2.2. For each state i// in a dense subset D ofH^in, o linear functional rjoii]/) in the algebraic dual of T) can be defined 
such that V0 e T> 

?7nWM lim ^ 

.^oJ^df(Qk''(M-)|Q)^,„ 

where Q e "TYa:,-,, is called a reference vector Therefore we can define a new inner product on the linear space of T]a(ip) via 
{r](>//)\T](4>)}n ■— ^nWW- The resultant Hilbert space is denoted by "Ha 

It turns out that the group averaging Hilbert space 'Ha is consistent with the physical Hilbert space defined by DID in a certain 
sense IfTSl : 

Theorem 2.1. We suppose zero is not a limit point in cr''''(M) and o"''''(M) = (which relates the argument that there is no state 
without physical interpretation). In addition, if we have any one of the following conditions 

1. there exists 6 > such that each (^Hm — l^m^^^) continuous on the closed interval [0, 6]. 

2. there exists 6 > such that each p'jf^ is continuous at A - Q and is differentiable on the open interval (0, 6). 

3. there exists 6 > such that N""^ is constant on the neighborhood [0, 6). 

Then there exists a dense domain D in 'Hkih, such that for some certain choices of reference vector Q the group averaging Hilbert 
space 'Ha is unitarily equivalent to the absolute continuous sector of physical Hilbert space 'H^'i^. 

See IfTSl for the proof of this theorem. Here we see that the reason of taking the limit e — > in Definition 12.21 is to make 
the desired connection between the group averaging Hilbert space 'Ha and the absolute continuous sector "H^^^^ in physical Hilbert 
space. For the pure point sector "H^'f,,, one should rather solve the eigen-equation Eq.( l2.12b in the kinematical Hilbert space 'Hkik- 
For the case of LQG or AQG, many eigenstates in 'HKin have been found, which correspond to the degenerated geometry, e.g. the 
spin-networks with valence less than 4. 

It is remarkable that all the physical models simpler than quantum gravity tested in ifTSi satisfy all the assumptions in Theorem 
12.11 It means that the group averaging technique in Definition 12.21 gives correct physical Hilbert space (the absolute continuous 
sector) for all those physical models. Therefore in applying the group averaging technique to the AQG master constraint operator 
M, we will assume that the master constraint operator M defined in Eq.( l2.11b also satisfies all the assumptions in Theorem 12. II in 
order to obtain the correct physical Hilbert space consistent with DID. 
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We will see in Section[3]that we can proceed the practical computation for the group averaging inner product 

r df <.Ak''^^"'Vk/« (2.13) 
Jb 

by using a standard skeletonization procedure ||231 . It turns out to be more convenient to use the coherent states in the skeletonization, 
which will be briefly reviewed in the next subsection. Finally, we will obtain a path-integral formula for the group averaging inner 
product. 



2.3 Coherent states in algebraic quantum gravity 

Before we derive the path-integral formula from the master constraint M, we brief review the definition of coherent state in LQG, 
in order to use the method of coherent state path-integral in the next section. Since we are considering an algebraic graph y with 
cubic topology, it is shown that the LQG coherent state is qualified as a proper semiclassical states in the kinematical Hilbert space 
'Hfcin- Its semiclassical expectation values for the geometrical operators and the master constraint operator reproduce GR correctly 
as ;i ^ lfT2ll9lfT3l. 

When we consider the semi-classical analysis of the above quantum framework, we should choose a embedding X of the algebraic 
cubic graph y into spatial manifold E, such that y is dual to the sub-2-complex >S in a cubic partition of S. Thus for each X(e) there 
is a face Se & S which intersects y transversely only at an interior point of both Se and X{e), for each x e 5^, we choose a path 
Pe{x) which starts in begin point of X(e) and along X{e) until pe and then runs within Se until x. We denote the path system for all 
X € S e hy "Ps ■ Given the data (y,X,S,'P), for a classical G-connection A and a Lie(G) valued vector density E of weight one, as a 
pair of conjugate variables, we define the holonomy and the gauge covariant flux for each edge e e E{y): 

hie) := !Pexp I A 

JX(e) 

Pie) := eabcdx" Adx'' hipAx))E'ix)hipAx))-\ 

And we define the dimensionless quantities pj - -^Pj - -^^Txir jP), tj - -/ (Pauli matrix)j for G=SU(2), and is a parameter 
with dimension [length]. By using the fundamental Poisson bracket between Aj, and E", we can check that the holonomies hie) and 
fluxes Pjie) form the following Poisson algebra lfT2ll : 

{hie), hie')} = 
[pj(e),hie')] = i-^6e,e'^hie') 

\pjie),pkie")\ = -i^Se,e'fjuPiie') 

The complexifier technique in LQG is motivated by the coherent state construction for compact Lie groups fTT]|. For separable 
Hilbert space, the coherent state transformation (which will be seen shortly) defined on a compact Lie group in 1111 is a unitary 
transformation from L^iG) to the BargmannSegalFock representation L^iG^) H HoliG^), where HoliG'^) is the space of holomor- 
phic functions on the complexified Lie group G^. We first proceed the construction on one copy of "He, which is isomorphic to the 
separable Hilbert space L?iG) with the Haar measure on the compact Lie group. We consider the phase space Me - T*G and define 
a, so called, complexifier as a positive function on the the phase space Me'. 

Ce f5''piie)pjie) (2.14) 
This complexifier complexifies the holonomies hie) by the following canonical transformation: 

8(e) = 2 {Ce,hie)\„^ = e-'^'^'^^^'^^e) 

n=l 

gie) belongs to the complexified group G'^ ^ T*G, which provides a complex polarization of the phase space Me- And (^gie),gie)^ 
are complex coordinates on the phase space Me- 
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Then we come to the quantization, we define a quantum complexifier Ce by 

Ce := ""ffPAe^pjie) = tAe (2.15) 

where pj{e) - iteXyi and are the right-invariant vector fields on G, A^, - d'-'X^X'j/A is the Laplacian on G, and tg - £^/al is 
a dimension-free classicality parameter The coherent state on edge e is defined by a heat kernel evolution of the delta function 
followed by an analytic contunuation: 

'^L.(Ke)) := le-'fSh'ieiiKe))] = [e'-^'^^^^w (/^(e))] 

= Yj^lje + 1) e-'-i'^'-^^'^I^Xh [gie)h(e)-') (2. 16) 

where g{e) is the complexified holonomy g{e) - e '^'*'^^'"j^^/i(e) e G*^. And the above heat kernel evolution followed by analytic con- 
tinuation h'{e) — > g{e) is called the coherent state transformation from representation to the BargmannSegalFock representation. 
Note that the above analytic continuation is trivially carried out, because the representation matrix elements of a compact Lie group 
G is the polynomials of h{e) in its fundamental representation. The label g{e) in the coherent state ^/r^'^^^ denotes the point of the 
complexified group G"' - where the coherent state is supposed to approximate. Moreover, it is easy to check that the coherent 
state (A^'j^) is an eigenstate of the annihilation operator g(e) with eigenvalue g(e) 

S^'^ -■-I.Mr ^^-'Hn, = S(e)^m = S(e)>f';,gy (2.17) 

It is important for us that the normalized coherent states iff^ form an over-complete basis in "H^, i.e. 



f 



g 

1 



dgie) l<Aj(,))<<Aj(,)l = 1, dg^ -dfinWd^p + oin (2.18) 



The coherent state on the fuU kinematical Hilbert space 'H^in is defined by a infinite tensor product 



^; = n ^;V) (2.19) 



eeE(y) 



The label g means that the coherent state i^^ is peaked at the phase space pointy = g {{h(e)]gf.E(y), ip iie)]eeE(y)) in the graph-dependent 
phase space Mx(y) - [r*G]'^*^" IT?! . For the operators Pol (^{h{e)]eeE{y), {pj{e)}eeE{y)j which is the polynomials of the holonomies 
and fluxes, the coherent state iff'g demonstrates the following semiclassical properties lfT2l : 

eeE(y), {p j(e)} eeE(y)) I lAj) = Pol {We)}eeE(y), {p j(e)} eeE(y)) + o(t) (2.20) 



More importantly, the semiclassical behavior of non-polynomial operators e.g. the volume operator V(R) and the master constraint 
operator M are also analyzed by the expectation value of the coherent state iff'^, by employing the semiclassical perturbation theory 
||9] [13] . The results for the volume operator and the master constraint operator on a cubic graph ar^ 

(.a; I y(/;) I (A;) = v(R)[g] + o(t) 

(iA;|M|«A;) = M[g]+o(t) (2.21) 

which demonstrate the correct semiclassical limits of both operators on a cubic graph ||9][l3|. This result shows that the coherent 
states on a cubic graph iff' is qualified to be the semiclassical states for AQG. 



o{t) is the quantum fluctuation which vanishes when ti 
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3 The path-integral of the master constraint 

Now we consider the practical computation for the group averaging inner product for defined in Definition 12.21 and derive a path- 
integral formula from the master constraint operator M on "Hkui- 



rcc 

J — CO 


dr 


^/ exp[/T(M- 


4f') 

J| /Kin 


/>CO 

J—oo 


dr 
2k 


|q |exp [;t(M - 


J 1 / Kin 



(3.1) 



where /,/' are kinematical states in a dense domain of 'Ka:,,, and M is the master constraint operator defined in Eq.( l2.11b . 
We multiply both the numerator and denominator of Eg. ( 13. II ) by the (infinite) constant 



J [DA(u)]6{^j'^ 



du A(u) — T 



(3.2) 



where [DA] is a measure on the infinite dimensional space of continuous functions {/l(M))«e[-r,r]- Note that here the integral 
du A(u) can be written as a Riemann sum lim/vi^oo jj Yjn=-N ^" ^^'-^ '^n = A(nT/N) when — > oo, provided A{t) is a contin- 
uous function. 

Insert this into Eg. ( 13. lb and perform the integral J dr 



lim 



LlT.(f\Mir{M-4f'),>„ 


J[DAiu)] 5(X-^j.dM A(t)-Tj 


/ [DA(u)] (f 


p[/T(M-e)] q) 
exp [ df A(u) (M 


J[DA(u)] J^J^j, du A(u) - rj 

-4f') 



/ [DA(u)] (q. |exp [ dM ^(m) (M - e)]| Q 
Therefore the central computation is for the matrix element 



(3.3) 



/ 



exp [ J du A(u) (M - e)] 



/' 



(3.4) 



by given a continuous function A(u) on the interval [-T, T]. By write the integral du A{u) as a Riemann sum: 



/ 



exp[ J duA(u){M-e)\ 



Kin 



exp[ ^ ^^(M-e)] 



Y] exp[^^(M-e)] 



Kin 



Kin 



Then we insert in the resolution of identity with coherent states Eg. ( 12. 18b on 'Ka:,„, 



(3.5) 



(3.6) 



we obtain 



Here dg„ = n«e£(y) d/i//(/z„)d^/9„/f^ is the Liouville measure on the graph phase space r*G'^*''". 

When we take the limit — > oo, we could approximate the matrix elements (as we usually do in the standard textbook-derivation 
of path integral formula) 



exp[|^„(M-e)]^;„)^^^^.(^;^ 



1 + ^^(M - e) 



'''^Kin 



(3.8) 
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Here A{t) is continuous function thus is bounded, so jjA„ is very small. Of course such an approximation still relies on non-trivial 
assumptions on the properties of the Master constraint operator (e.g. some assumption described in ||27l ). 

First let's compute the single-step ampHtude (lAgJl + '-^^ni^ ~ Here the overlap function {>jy\^\4''g^ ^)Km is sharply 

peaked at gi - in a Gaussian fashion (with width yft). Thus in the semiclassical limit f — > we have 



1 + ^^(M - e) 



'Kin 



iT <iA' IM-elfA' >^ 



1 + —A„ 

N 



lim 



iT 



ikn i + -^(mm-6) <^;j^;,_,). 



In the second step, we use the fact that as f ^ 0, the overlap function (i^^ Ji^^^^ )jf„, is proportional to a delta function (5g,(§,_i). In the 
last step we use the fact that the master constraint operator has correct semiclassical limit Eq.( |2.21t . Therefore, we have seen that 



1 + ^^(M - e) 



gi-\ 



Kin 



\ + '^A„{M[gi\-e + tF'(gi,gi-i)) 



N 



{4''g}l('ljKin 



where F'(gi, is the fluctuation of the master constraint operator M with respect to the coherent states ifr'^. and i^^^^ . Therefore as 
we formally take that limit ^ oo, we have approximately 



- exp 



'■^A„ {M[gi]-e + tF'(gi,gi^i)) 



N 



Then let's evaluate the overlap function (i/;-^J0'^, _);f,„ = nee£(y)<'/'j5,(e)l'A'^,,_,(f))««/l|i/'|;,(<,)ll H'/'^.^itf)!! where 

CO 

^'g(e) We)) = ^^j' + l)e-'j'^j'^'''\j^{g(e)h-\e)) 

2>,=0 

is the complexifier coherent state on the edge e. If we set n = 2^ -H 1 € No, on one edge 



(K>\K<-,>Kin = 



2sinh(z,;,_i) 



n e e '•' 



(3.9) 



(3.10) 



(3.11) 



where cosh(z,y) = tr(g]gj)/2. For the norms of i^^ ,, the same formula holds just with replaced by z; such that cosh(z,) - tT(g .gi)/2. 
Therefore, 



<<i''g.\4''gjKi„ = n 

eeE(y 

" n 



2sinh(z,-,_i) 



.in-Z'i nzij- 



VI sinh(z,) sinh(z,_i)| 



gtn^HgnZi-i 



ZiJ-\ 



eeEiy) 



sinh(z/,,_i) 



^j\Zie''l'\^j\Zi-xe^^l'\ 



pr VI sinh(z;) sinh(z;_i)T ^,,2_^_ 1 ,2_ 1 .2 ^y, 
elEiy) VkiZ/-il sinh(z/,,-_ I ) 

where we have neglect a term of order o(f°°) in the second step, and we have used the Poisson resummation formula 

Yfins)^-Y r dae^^'^l^fix) 



(3.12) 



(3.13) 



with s - V^c- 

On the other hand, we have the relation: 



cosh(z,) = tr(g,'g;'")/2 = \x{h\e-'^P'-'l^e-'P'-'l^h-:)ll 
= tr(e-'>'-^)/2 = cosh(;7,) 



(3.14) 



which gives |z,| = /?,■ = ^|pr~Pi■ 
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Then we can apply the above result for the overlap function to compute the matrix element (as — > oo) 



/ 



exp [ J du A(u) (M - e)] 



/' 



Kin 



N 



- lim 



n 

k=-N+l 



Inn 

^ eeE(y) i=-N 
Zk,k-l 

sinh(zt,t_i) 



dhjtPpi sinh(p,) 



-Pl/t 



e~k.k-\l' exp 



'^A, {M{gk\-e + tF\gk,gk-x)) 



N 



f(gN)f'ig-N) 



(3.15) 



where /(g) := {^'^ \ f}Ki„. 

Furthermore, we define the Lie algebra variables 0^ such that hi^ = e"' Zk,k-i can be computed by 



cosh(zk,k-i) 



-mk-ihle-'^P'^P^-^^-^l^) 



= coshj^ 



iPk+Pk-l) .(dk- dk-l)-\2 
h I ' 



Therefore, 



2 _ 2 _ 1 2 



~[(Pk - Pk-xf + {Ok - Ok-xf + 2i(pk + Pk-i) ■ {Ok - Ok-x)] 



(3.16) 



(3.17) 



Then we insert the result back into Eq.( |3.15l l and obtain 

I-T 



exp [ ^ du A{u) (M - e)] 

fc/ n fi ^ 



f 



Kin 



sinhipi) 



(Pl«+P«)/2' 



cE£(y) i=-N 



n 



j(:=-A'+l 



sinh(zt,i_i) 



.figN).f'{g-N) 



X exp < -i 



Xpk +Pk-x) 
2t 



■ (Ok - Ok-x) - UiPk - Pk-xf + (Ok - Ok-xf] + '^Ak [Mfo] - e + tF'(gk, gk-x)] \ (3.18) 



In order to contact with the path-integral formula with a classical action on the exponential, we make the following approximations 
forEq. dTTSlR : 

• We assume the fluctuation F' of the master constraint operator M with respect to the coherent states ij/'g is small and negligible. 
It is a non-trivial assumption for the property of the master constraint operator M, i.e. we should design a certain operator 
ordering in the definition of the self-adjoint master constraint operator M, such that the fluctuation F' is small and negligible. 
For simple system Hke the ordinary free quantum field theory, such an operator ordering is nothing but the normal ordering of 
the creation and annihilation operator, which results in F' = 0. 

• We only count the paths such that the second order terms (Ap)^ and (AO)^ on the exponential is negligible. Note that the 
overlap function (0''i;J'Ag,_,)A:m sharply peaked at the point gi = g/_i with width sfi. Thus (Ap)^, {Aff}^ ~ t and have to be 
neglected if we want to have a classical action on the exponential. It is consistent if we want to have a classical kinetic term 
containing the time-derivative of 9 after the continuum limit. 



'These approximations essentially remove all the /"^''-order contributions on the exponential, and keep only the t '-order contributions. 
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With the above approximations. Eg. ( 13. 181 ) is simplified to be 



exp[ J duA(u){M-e)\ 

= if J n fi ^ 



f 



Kin 



sinh(/:>,) 



-(piN+PN)nt 



n exp^-i 

k=-N+\ 



eeE(y) i=-N 

(pk +Pk-\) 



n 



Zk,k-\ 



.figN).f'{g-N) 



iT 



It 



■(ek-0k-i) + —Ak [M[gk]-e] 



= j [Dei(x")Dpp^)] exp I - J dx^ 2 p'^jd.oei + i dx" Aix") [M[pp\ - e] 



eeEiy) 



f[p''jiT),9i(T)]f'[p%-T),ei(-T)] 
where we have relabel the coordinate u - and the functional measure is formally defined by 



[dO^AADpP^ := hm f] f] 



dhid^Pi sinh(/?,) 



N^ca 1111 p 

eeE(y) i=-N 



n 

k=-N+\ 



Zk,k-\ 



sinh(zt,n) 



(3.19) 



(3.20) 



Then we insert Eg. (13.191 1 back to Eg. (13. 3b and perform the functional integral of DA (we consider the numerator) 
lim ^ {DA{u)] {f exp [ ^ dt A(u) (M - e)] /'| 

= lini J [DA(x^^)] J [Dei(x°)Dp'j(x°)\ expj-^ J dx° P'jd.^ei + i J dx° A(x°)[M[p'j(x\ei{x^^)] - e^^ 

f [p^), ei(T)]f' [p%-n 0i(-T)] 
= lim f[Dei(x°)Dpp'')] expl-- f d/' J] p)d_,oei\ 6 {M[pp\ 0i(x°)] - e) 



f[p''j(T),0i(T)]f' [ppT),6i(-T)] 



(3.21) 



Now we consider the delta function 5^M[p^(x''),0^(jic°)] - ej for each jc". If the graph y here was a finite graph with a finite 
number of vertices, thus the sum in master constraint 2i e£(y) would be a finite sum. However, even we consider an infinite graph 
with infinite vertices, i.e. the master constraint has a shape M = C/C/ with some constraints C/, where the sum in the expression 
is a infinite sum and J is a infinite index set. We can truncate the sum and define a partial master constraint Mw - Tjiew CiCi where 
W c J is a finite set. Then it turns out that the group averaging inner product by using the partial master constraint Mw converges to 
the group averaging inner product by using full master constraint M as W — > X under some certain technigue assumptions (see lITSl 
for the details). It means that for a infinite graph y we could define a delta function for the partial master constraint. 



MM-e):=5(M„,-e) such that lim 6w(M - e) ^ 6(M - e) 



(3.22) 



in the sense of distribution. 

So here we only consider the following type of integral with a finite integer and any test function / continuous in a neighbor- 
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hood of the constraint surface: 



lim I dxidx2 • • • dxN S( } - e) /(:?) 
= hrn I dx2 • • • dx^ , f(xi = Je- S^j-^'-^') + /(-^i = ~ ~ ^i^2-^'-^') 

= {1™/ 



/(;?=o) 



= hinA^(e)/(x = 0) 



Mm Nie) I dxidx2 • • • dxf^ 6{Xi) /(^ 



where ATCe) diverges in the limit e ^ 0. 

So far for the integral on the space of ^, if we parametrize by the variables {^yl^i (M > N) and consider the following 
integral: 



lim I dri---d;M^(y Xi(^>i(^l-e)/(^"3• 



Then we make the changing of variables 



tk for k€{l,--- ,N} 

tk tk for k € {N + I,- ■■ ,M} 

and denote the Jacobian to be det{dx/dt) and assume it is nonzero, continuous and bounded around the constraint surface. Therefore 

Mm I d^i • • • dtM Xiiijxii^ - e) f{h, ■■■ Jm) 

= hm I dxi • • • dxArd?A,+i • • • dtM , ^(Y, XjXi - e) /(xi, • • • , X;v, ^w+i, • • • , ^m) 

e-»o J detidx/ot) j-^ 

= limyV(e) I dxi • • • dxt^dt^+i ■ ■ ■ dt/n , ,„. O <5(xi) /(xi, ■■■ ,xn, tN+i, ■■■ Jm) 
e->o J detidx/ot) 11 

= lim Nie) J dfi • • • dfM Y\ ^(^t) fiti Jm). 



i=l 

Therefore we obtain the result: 

N 



lim 6i y x,(?lx,(?) - e) = lim Nie) Ff <5(x,) (3 .23) 

e-»0 ' e— >0 i A 



(-1 1=1 



Since the master constraint in our case has the following expression: 

Gj^V Gj^V Dj^l, Djy Hy Hy 



where Gy, and H here are the original Gauss constraint, diffeomorphism constraint and Hamiltonian constraint in their regularized 
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forms, then we obtain the following path-integral formula 



/' 



hm J [DA{u)] |/ exp [ J dt A(u) (M - e)] 
limyV(f) r [£>0i(x°)Dp5(x°)] exp -| f dx" ^ /,}5,o0> \ f[p'^j{T),0i(T)]f' [p''ji-T),9ii-T)] 

limyV(e) J [z)0i(xO)Z)p5(xO)DA;:(x")DA^/(x'')DA^(x°)] / [p}(r), 0^(r)]/' [pj.(-r), 
exp|-J r dx" 2 /?55,o^ + 2 {AiGj,, + NiDj,.+NM.) \ 

[ Lee£(y) veVir) J j 



(3.25) 



where the factor limf^o A/^(e) can be canceled out from both denominator and numerator of Eg. (13.31 1. Here the path-integral measure 
is defined formally by 



[Ddi(x^)Dp'j{x")DA{,(ADNi{x^*)DN(x^) 



:— lim 



n n 

1=-^ eeE{y) 



dhi(Ppi sinh(;?,) 



-(P-A.+P«)/2' 



n 

k=-N+l 



Zk,k-l 



sinh(zi,n) 



n=-N+\ veV(y) 



In the end we write down this path-integral representation of the group averaging physical inner product 

<'7(/)l'7(/')>n = 



Zy(/,/') 



Zy(D, Q) 



(3.26) 



with 



ZyifJ') = J [De^Ax'')Dp)U')DMix'')DNi{x'')DN{x'')\f[p){T\^^^ 

expJ-^ r dx» ^ ;?}5,o0/+ ^ (a/:G,>+A?/d,> + A?,//,) I. 



(3.27) 



There are a few remarks for the path-integral formula Eg. (13. 25b : Motivated by the following relation between ^(e) - e 'P'*^'^^'^'l^h{e) e 
j*Q\E(y)\ ^jjjj j}jg classical phase space variable (A^, E") € M via the embedding x^ : y ^ 2 



/z(e) 



]pj{e) = -tr(^ ^ dx" A dx*e„i,/i(p,(x))£^(x)/z(p,(x))-') 



Then obviously Eg.( l3.25l l is in analogy with a path-integral formula on the spacetime manifold M when we take a certain continuum 
Umit (recall that / = fpla^): 



J [DAiDE"pAjDN"DN] f[Ai,E''j]^f'[AlE'^]_^ 
- ^ J d-^x {E"jdMi - J^^Gj - N"Ha - NH) 



X exp 



(3.28) 



It is obvious that the integral on the exponential is nothing but the canonical action of GR in terms of su(2) connection and electric 
field variables. And the path integral measure is the continuum version of 



[D9i(x^)Dp''j(x^^)DAi(x°)DNi(x'-^)DN(x^)] 



:= lim 

N-*oo 



n n 

i=-N eeE{y) 



dhid^Pi sinhipi) ^_(„2 , „2 



(/'-A.+/'«)/2( 



n 

k=-N+\ 



Zk,k-l 



sinh(zi,n) 



N 1 r 

Y] Y] d^A,.,„d^N,,,m,„ Yl n ^''^^P"^ 

n=-N+l i-eV(y) \[n=-N+l veV(y) 
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Thus the naive Ansatz of the path-integral formula Eg. (1 1.4b is reproduced by Eg. (13. 28b up to the measure factor The spatial volume 

7 II 

local measure factor Vs appeared in the path-integral measure is because the master constraint is of the type M = 2/ KjjCiCj 
(recall Eg. ( 13. 24b ). with a phase space dependent Kij. However for example, if we define the master constraint by 

M = ^ Gj,,Gj,, + Dj,,Dj,, + HMv (3.29) 

instead of Eg. (13. 24b . the spatial volume factor v]^^ disappears in the path-integral formula. However, the definition Eg. ( 13. 24b is 
preferred because it has diffeomorphism invariant continuum limit. 

4 Conclusion and discussion 

In the previous analysis, we employ the group averaging technigue for the AQG master constraint M on a cubic algebraic graph y, 
and derive the expected path-integral formula for the group averaging physical inner product Eg. ( 13. 26b . which is a certain analog 
of the Feynman-Kac formula in guantum mechanics |[26l . The derivation uses the semiclassical tools developed in |[9][l2|, thus the 
final path-integral formula is in the form of a coherent state path-integral. Although the resulting path-integral formula relies on 
non-trivial assumptions and approximations, our analysis suggests the existence of a path-integral formulation directly relating the 
master constraint guantization for GR. 

Another implication of the previous analysis is that, we suggest a more direct path linking canonical LQG with a path integral 
guantization of GR, comparing with the approaches we suggested in Il24l . The procedure in Il24l is the following: We start with the 
canonical theory of GR (the Hoist and the Plebanski-Holst formulation) on the continuum, and perform the formal guantization in 
its reduced phase space. Then we derive from the reduced phase space guantization a formal path-integral of GR on the continuum, 
which is formally proved to represent the physical inner product of GR. 

However the derivation in the present paper start from Dirac guantization in terms of a concrete operator formalism. We start 
with the guantization of GR via AQG on an infinite algebraic cubic graph 70 and its self-adjoint master constraint operator M. The 
physical inner product is formulated in terms of the group averaging inner product of the self-adjoint master constraint M. The 
skeletonization procedure for the group averaging of M arrives at a path integral representation Eg. (13. 261 ) of the group averaging 
physical inner product from the operator formalism. Interestingly the skeletonization procedure in the derivation shows some certain 
similar structures between the group averaging physical inner product and the spin-foam guantization (on a hypercubic spacetime 
lattice). Thus it could be a future research about derive a spin-foam vertex amplitude from e.g. Eg. ( 13. 151 ). 

The resulting path-integral formula in the present analysis is in a canonical form Eg. ( 13. 25b . i.e. the action on the exponential 
is the canonical action of GR and the path-integral measure d/i,d^/7, is the Liouville measure on the phase space 7"*G'^^^^'. It is 
interesting to consider deriving from this canonical path-integral to a covariant path-integral with original covariant field variables, 
where we may have to make a certain unfolding of the path-integral measure Il24l and use the Henneaux-Slavnov trick |[29| . From 
this procedure it would be interesting to see if the SO(4) or SL(2,C) holonomies of the spacetime covariant connection fields can 
appear in the final formulation. All current spin-foam models 1201 start from a discretized path-integral formula in terms of the 
SO(4) or SL(2,C) holonomies and the Haai" measure of SO(4) or SL(2,C), but the result here Eg. (l3.25l ) only display a SU(2) Haar 
measure d/z(e) product with the Lebesgue measure dAy standing for the spacetime connection field. 
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'"We can also start from LQG with a chosen graph y, the deiivation for the path-integral formula is the same as we did in Section[3] 
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